Open End Correction for Flanged Circular Tube 
from the Diffusion Process 

Naohisa Ogawa 

E-mail: ogawanao@liit.ac.jp 

Hokkaido Institute of Technology, 7-15 Maeda, Teine Sapporo 006-8585 Japan. 
Fumitoshi Kaneko 

E-mail: toslii@cliem. sci . osaka-u. ac . jp 

Department of Macromolecular Science, Graduate School of Science, Osaka 
university, Toyonaka Osaka 560-0043 Japan. 

Abstract. In the physics lessons on waves and resonance phenomena in high school 
and college, we usually consider sound waves in a tube with open or closed ends [T]. 
However, it is well known that we need a tube with open end correction AL. The 
correction for a flanged circular tube was first given by Rayleigh [2j and experimentally 
checked by several authors j3]. In this paper, we show the different methods of 
obtaining the end correction for a circular tube by a diffusion process. 
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Figure 1. Moving piston in a tube connected to half infinite space filled with 
noncompressive fluid. 

1. Introduction 

A standing sound wave in a tube with an open end has an antinode bit outside the end 
point. Such a small distance is called the open end correction (hereafter abbreviated 
as OEC). The end correction AL was calculated for an infinite flanged circular tube by 
Rayleigh [2J, and after that, H. Levine and J. Schwinger calculated the unflanged end 
correction as a function of wavelength 0] . OEC for an infinite flanged circular tube as a 
function of wavelength was determined by Y. Nomura et al. by the radiation impedance 
method [5]. The recent application of OEC was given by M. S. Howe [6 J (and related 
references therein) . 

The idea of obtaining OEC as discussed by Rayleigh comes from the fluid 
mechanical consideration [2], [6]. Let us suppose a half infinite space with a connected 
circular tube filled with nonrotational and noncompressive fluid, and assume a piston 
moving at a speed of V in the tube as seen in fig. 1. Then, the total kinetic energy of 
the fluid is given by the definition of OEC. 

l -pj(y<P) 2 o?x= l -pV 2 S(L + AL), (1) 

where p is the fluid mass density, is the velocity potential, S is the cross-sectional 
area of the tube, and AL is the OEC. The integration region is taken toward the piston 
inner circular tube and half infinite 3D space. The region with a nonzero velocity is not 
only in the tube but is also spread around a branching bay with a speed lower than V. 
Such a region contributes to OEC. 

In this study, we calculate OEC quite differently. As shown in fig. 2, we place a 
circular well at x — y — with a radius a and a depth L, and consider the particle 
diffusion of density n(x, y, z). We put a particle bath with a density n s at the bottom of 
the well; we thus have n = at infinity. By using the above boundary conditions, 
we approximately solve the static diffusion equation under the condition that one- 
dimensional diffusion is occurring in the well. By using a solution for a diffusion field, 
we define OEC and show that its value is 8a/37r, consistent with that given by Rayleigh 
and Nomura et al., when the tube radius is small enough compared with the wavelength 
a/A « 1 |5]. 
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2. Estimation of diffusion in a well 



In the well, we approximate the diffusion process by using a one-dimensional diffusion 
equation (see fig. 2). This is very simple, but is essential for our further discussion. At 
the top of the well, the space is widely opening and it is therefore possible to suppose 
that n top = n(z = 0, r < a) = 0. Then, the constant diffusion flow comes up from the 
bottom. 

J 

n{z)=n s --{z + L), J = — ^. (2) 
Our purpose is to correct this flow formula. 

We take nt op ^ and then obtain the diffusion field and constant diffusion flow. 

n{z) = n s - J -{z + L\ j> = !h-^. (3) 

na 2 J' is the number of particles supplied per unit time from the well to the open region 
while leaving n top undetermined. Then, we solve the diffusion equation in open space 
with the boundary conditions 

n(r, z = oo) = n(r = oo, z > 0) = 0, (4) 

dn _ J • • • r > a, z — , . 

dz = \ -J'/D ■■■ r<a,z = 0. ^ ' 

The additional condition 

n(z = 0,r < a) = n top (6) 

is given below using the solution of the diffusion equation in an open region. Then, two 
diffusion fields in the inner and outer wells are connected, and the diffusion equation is 
completely solved. 

3. Diffusion in open region 

We consider the diffusion equation An = in the open region z > 0. The boundary 
conditions are flU) and (JSJ). 

By supposing an axial symmetry, we obtain the following static diffusion equation 
for the diffusion field n: 

d 2 n 1 dn d 2 n 
dr 2 r dr dz 2 
Then, we have a solution in the form 

roo 

n{r,z)= d\f x J (\r)e- Xz . (8) 
Jo 
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Figure 2. Circular well with depth L. The particle server is at the bottom. Particles 
diffuse from the bottom of the well into the outer region. 



On the plane z = 0, the diffusion flow in the +z direction is given by 

j( z = 0) = _ Ddzn \ z=0=D \f x J (\r)dX. (9) 

J 

fx should be selected to obtain the boundary condition (jSJ). 
The following formula solves the problem. 

WW ;;;<;: m 

From (j2J), and ( fl~0l) . we determine the function /a given as 

A = ^W (ii) 

Then, we obtain the following solution for n: 

n(r, z) = ^- T ^ J x (aA) J (Ar)e- Az . (12) 

Next, we must determine n top = n(r < a, z = 0). From (j!2p . we obtain 

n(r,0) = ^ r -.h(x)U-x) = ^N(r/a). (13) 
D Jo x a D 

The function N(r/a) is expressed by the following hypergeometric function: 

[ F(l/2,-l/2,l;(r/a) 2 ) ••• r<a 

N(r/a) = < 2/tt • •• r = a (14) 

( F(l/2, 1/2,2; (a/r) 2 )/(2r/a) ■•• r > a. 

The form of this function is given in fig. 3. The smoothed curve is the function N(r/a). 
To obtain a consistent solution with a diffusion field in the well, the matching condition 
(JO]) is required at (r < a, z = 0). However, as shown in figure 3, the distribution at 
r/a < 1 is not constant. This shows that our solution ffT2|) is not consistent with the 
inner-well solution, since our approximated solution (j3J) counts the diffusion out in the 
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Figure 3. Density distribution on surface z = 0. At r > a, we have a nonzero density. 
This shows that particles diffused from the well and distributed around it. The square 
line shows the mean N[r/a) for r < a. 



r-direction. To further achieve what in this approximation, we approximate the form of 
N(r/a) by its mean in the region r/a < 1. 

<N>=% T drrF(l/2, -1/2, 1; (r/a) 2 ) = K = ^- ~ 0.8488. (15) 
a 2 Jo 3n 

This mean is shown by a solid square line in fig. 3. Then, an additional boundary 
condition (jHJ) leads to 

ntop = -jy J '- ( 16 ) 

From <$\) and (JIB]) , we obtain 

J ' = (L + aK)/D' (17) 

This result should be compared with the second equation of (j2J); it shows that the 
effective depth of the well is slightly larger than L and that OEC is defined by 

AL = aK = — . (18) 
3n 

This is the same OEC for sound in a long- wavelength limit [5]. The point z = AL 
upwards of the well is the place where the density n vanishes effectively. 
Note that, from eq. (JIB]) OEC is defined by 

AL = (19) 
J' jD ' V ; 



where 



ira 2 Jr<a 



< n >= — - / n(r, z — 0) dS. 
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4. Discussion 

Why do we have OEC for sound in a long-wavelength limit by our method? We discuss 
this point and show that the definition of OEC by the radiation impedance method is 
quite similar to that given in this paper. 

First we review the radiation impedance method of obtaining OEC. Let us put a 
disk as a sound source oscillating in the z direction and having a radius a on the surface 
z — 0. The sound emitted from the source satisfies the wave equation, 

(A " 7*W )p = °' (20) 
where p is the pressure field and c is the speed of sound. 

The linearized Euler equation is necessary to consider the speed of sound source: 

'f - < 21 > 

where p is the mass density. 

From the solution of eq. (!20|) and ( 12TI) . we obtain the fluid velocity u. 
The boundary conditions for p at the surface z = are as follows: 

dp = | (u z = 0) ■ • • r > a, z = 0, 

dz \ -iu p u z (source) • • • r < a, z = 0, 

where the time dependence is supposed to be exp iujt. 

From the emitted sound p(x,t), we obtain the force of reaction F to the sound 
source as 



F= / P{z=0) dS. (23) 



r<a 



Then, the radiation impedance is defined by 



F - -ipu Ir <f^ dS . (24) 



Uz (source) d z p( r<a ) 
From the parallelism to the mechanical impedance, the additional inertial mass due to 
the radiation is defined using the maginary part of radiation impedance. 

5m = -9f(Z). (25) 
to 

This additional mass is supplied by OEC as 

5m = p AL (vra 2 ). (26) 
So, we obtain OEC by using the formula 

AL = m. = -S(^), (27) 
pixa 2 uj o z p {r<a) 



where 



<P>= — ^ / P( z =o)dS. 

7Ta z Jr<a 
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Now, let us check the similarity of this OEC in a diffusion process. The long- 
wavelength limit is given as 

A — > oo, oj — > 0. 

Then eq. ( 120]) leads to the Laplace equation. The linearized Euler equation ( I2TT) 
should be compared with the equation of a diffusion flow. 

J/D = -Vn. 

Therefore, we have the correspondences; n ~ p and J/D ~ icupu. The boundary 
condition ( 1221) corresponds to 

dn = { ■■■ r>a, z = 0, 

dz \ -J/D ■ ■ ■ r < a, z = 0, 1 ' 

which is the same as ([5]). 

The definition of OEC as eq. ( 1271) corresponds to 

Ai = = (29) 

d z n (r<a) J/D 

This is the same as the definition of OEC for the diffusion equation ( TT9l) . Therefore, 
the equations, boundary conditions, and definitions of OEC are the same in the long- 
wavelength limit of these two methods. This is the reason why we can calculate OEC 
of sound waves simply using a diffusion equation. The difficulty in obtaining OEC is 
quite simplified in the long-wavelength limit by considering the diffusion equation. 
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